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I. INTRODUCTION 

As hybrids of form factorsjtarton distribution functions and distribution amplitudes, the generalized parton dis- 
tributions (GPDs) Fp, 0, H, 0, @, 0| provide the most complete and unified description of hadronic structure (for recent 
reviews, see 0, 0, Ej] )■ Parton distributions parameterize the longitudinal momentum distributions (in an infinite 
momentum frame) of partons in the nucleon, while the Fourier transforms of form factors in impact parameter space 
describe the transverse coordinate distributions of the nucleon's constituents [Iol.[ll|. GPDs, on the other hand, simul- 
taneously encapsulate both the longitudinal momentum and transverse coordinate distributions, and hence provide a 
much more comprehensive, three-dimensional snapshot of the substructure of the nucleon. 

At the same time, GPDs have received considerable attention in recent years in connection with the so-called "proton 
spin crisis" (l2l. H3. H3| . Namely, certain low moments of GPDs can be related to the total angular momentum carried 
by quarks and gluons (or generically, partons) in the nucleon 0|. Combined with measurements of the quark helicity 
from inclusive deep inelastic scattering, knowledge of GPDs can thus unravel the orbital angular momentum carried 
by partons, on which little or no information is currently available. 

Typically, GPDs can be measured in hard exclusive lepton-production processes, such as deeply virtual Compton 
scattering (DVCS), 

e-(k)N( Pl ) — e-(k')N{p 2 ) 7 (q 2 ) . (1) 

Here, an electron (or muon) scatters off a nucleon via the exchange (in leading order QED) of a space-like photon with 
virtuality qf — (k — k') < 0, producing an intact nucleon (with altered momentum) and a real photon in the final 
state. (The four-momenta of the particles are denoted in the parentheses.) At the quark level, in the leading-twist 
approximation, the electromagnetic current couples to different quark species with strength proportional to the squares 
of the quark charges, selecting specific linear combinations of GPDs. Flavor-specific GPDs can be reconstructed by 
considering DVCS from different hadrons (protons and neutrons, for instance), and using isospin symmetry to relate 
GPDs in the proton to those in the neutron. 

On the other hand, different combinations of quark flavors can be accessed by utilizing the weak current, which 
couples to quarks with strengths proportional to the quark weak charges. In analogy with DVCS, these can be studied 
in neutrino-induced virtual Compton scattering, 

v (k) N (pi) — e-(k')N'( P2 ) 7 (q 2 ) (2) 

for the charged current, and 

v{k)N{p x ) — iy(k')N(p 2 ) 1 (q 2 ) (3) 

for neutral current reactions (and similarly for antineutrinos) . In particular, neutrino- induced DVCS can be more 
sensitive to the d quark content of the proton, in contrast to electromagnetic probes which, because of the quark 
charges, are sensitive mostly to the u quark. 

Because of the V — A nature of the weak interactions, one can probe C-odd combinations of GPDs as well as C-even 
(where C is the charge conjugation operator), and thus measure independently both the valence and sea content of 
GPDs. This has a particularly novel application in the case of the polarized distributions. The usual way to obtain 
information on the spin structure of the nucleon is through inclusive deep inelastic scattering of a polarized lepton 
from a polarized target. To separate the valence and sea spin contributions one could imagine measuring the C-odd 
polarized structure function by scattering neutrinos from a polarized target, which would be rather prohibitive using 
existing technology given the large volume of target needed to polarize. Because it is sensitive to both spin-averaged 
and spin-dependent GPDs, neutrino-induced DVCS could allow one to extract the spin-dependent valence and sea 
quark distributions using only an unpolarized target. 

The weak current also allows one to study flavor nondiagonal GPDs, such as those associated with the neutron-to- 
proton transitions in charged current reactions in Eq. <[2j) ■ The use of weak currents can thus provide an important tool 
to complement the study of GPDs in the more familiar electron-induced DVCS or deeply exclusive meson production 
processes. 

Recently, neutrino scattering off nucleons for neutral currents was discussed in Ref. [l5| , where the authors presented 
the leading twist behavior of the cross section for the deeply virtual neutrino scattering process. In Ref. [IB], the 
analysis was extended to charged currents, and the leading twist amplitude computed. The neutrino-induced hard 
exclusive production of D s mesons was considered within the GPD formalism in Ref. [ItJ • 

In this paper, we present a comprehensive analysis of the weak deeply virtual Compton scattering (wDVCS) 
processes in Eqs. ^ and |[3|), and give a detailed account of the charged and neutral current amplitudes and cross 
sections in the kinematics relevant to future high-intensity neutrino experiments [l8| . Some of the formal results in 
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this paper have appeared in a preliminary report in Ref. [19]. In Section HH we provide a detailed derivation of the 
leading twist weak neutral and charged current amplitudes using the nonlocal light-cone operator product expansion. 
We introduce an appropriate set of GPDs, which parameterize the wDVCS reactions. Together with the standard 
electromagnetic DVCS process, the latter are analyzed in Section ITTll where we also discuss the relevant cross sections 
and kinematics. Here we closely follow the analysis of Ref. [30], however, we only keep contributions up to the twist-2 
accuracy. Using a simple model for nucleon GPDs, which includes only the valence quark contribution, we estimate 
the cross sections for weak DVCS processes, and further compare the respective rates in neutrino scattering with 
those in the standard electromagnetic DVCS process. Finally, we draw some conclusions and discuss future prospects 
in Section [iVl 

II. WEAK VIRTUAL COMPTON SCATTERING AMPLITUDE 

In this section, we present a detailed derivation of the amplitudes for weak virtual Compton scattering. We begin 
with an analysis of some general aspects of the amplitudes, before turning to the specific cases of the weak neutral 
and charged currents. 

A. Generalities 

In analogy with the photon-induced DVCS amplitude, the weak virtual Compton scattering amplitude can be 
obtained by replacing the incoming virtual photon with the weak boson B, 

B(q 1 )N( Pl ) — ► 7 (Pa) , (4) 

where B = Z° or W . In the case of the charged W bosons, the initial and final nucleons will be different. In the 
Bjorken regime, where the virtuality of the initial boson and the total center of mass energy squared of the virtual weak 
boson-nucleon system are sufficiently large, namely — q\ and (p\ + q±) — > oo with the ratio xb = —Qi/ [2 (pi • qi)] 
finite, the relevant amplitude is dominated by light-like distances. The dominant light-cone singularities, which 
generate the leading power contributions in l/|?i| to the amplitude, are represented by the so-called "handbag 
diagrams" in Fig. 1. At leading twist and to the lowest order in a 8 , there are two diagrams which contribute, in which 
the (hard) quark propagator is convoluted with the (soft) four-point function parameterized in terms of GPDs. In 
addition, keeping the momentum transfer squared to the nucleon, t = (p\ — P2) 2 , as small as possible, one arrives at 
the relevant kinematics required to study DVCS. One of the methods to study the amplitude in these kinematics is 
based on the nonlocal light-cone expansion of the product of currents in QCD "string" operators in coordinate space 
[2p| . which we will employ in the present work. 




Figure 1: Weak deeply virtual Compton scattering of a weak vector boson B from a nucleon N, producing a real photon 7 and 
recoil baryon N' in the final state. In the leading "handbag" approximation (two rightmost diagrams), the boson scatters from 
a single quark in the nucleon, which emits the real photon. 

In the most general nonforward case, the virtual Compton scattering amplitude is given by a Fourier transform of 
the correlation function of two electroweak currents. In particular, for the standard virtual Compton process on the 
nucleon, both currents ( Jem an d ^em) are electromagnetic and the amplitude can be written as: 

n v M = if*xj d%e- i ^ +i ^y{N(p 2 ,s 2 )\T{J^ M (y)J EM (x)}\N(p 1 ,s 1 )) . (5) 
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Here p\ and p 2 are the four- momenta of the initial and final nucleons, and Si and s 2 are their spins. Similarly, for the 
weak process, with an incoming Z° or W boson and outgoing photon, the amplitude is: 

T% = if d 4 x f d 4 ye- i ^ x+i ^(N'( P2 ,s 2 )\T{J^ M (y)J^(x)}\N( Pl ,s 1 )) . (6) 



Here (x) corresponds to either the weak neutral current J^n or the weak charged current Jyyc Note 
that in the electromagnetic and weak neutral cases both the incoming and outgoing nucleons are the same, N = N' , 
whereas for the charged current case they are different, N ^ N' . 

It will be convenient in the analysis to use symmetric coordinates, defined by introducing center and relative 
coordinates of the points x and y, X = (x + y) /2 and z = y — x. The weak virtual Compton scattering amplitude 
then takes the form: 

T% = if d A X f d*z e -«(tt-fe)-*-Hte+*Wa <jv/ (p 2 ,s 2 )\T {J» M (X + z/2) (X - z/2)}\N (p 1)Sl )) . (7) 



Furthermore, in order to treat the initial and final hadrons in a symmetric manner, we introduce as independent 
momentum variables the averages of the boson and hadron momenta, q = (qi + q 2 ) /2 and p = (jp\ + p 2 ) /2, and the 
overall momentum transfer, r = pi — p 2 = q 2 — q\. Accordingly, we have: 

q 2 = q\/2 - t/4 with t = r 2 . (8) 

From the on-mass-shell conditions, p\ = M 2 and p\ = Mf , one then has: 

P 2 = \ {M\ + M 2 - t/2) and p ■ r = i (M 2 - M 2 ) , (9) 

where Mi and M 2 denote the masses of the initial and final nucleons, respectively. In the following we shall neglect the 
mass difference between the proton and neutron, and set M\ = M 2 = M(~ 0.94 GeV), in which case these relations 
simplify to: 

p 2 = M 2 - and p ■ r = . (10) 

After translating to the center coordinates in Eq. ([7]), (p 2 \ J M (X) \pi) = (p 2 \ J p (0) \p\) e~ l ( pi ~ P2 ^' x , and integrating 
over X, one can write the weak amplitude as: 

T% = (2tt) 4 (p x + qi -p2- q 2 ) 7&" , (11) 

where the reduced weak virtual Compton scattering amplitude is: 

Iff = ij d 4 ze^ z (N'(p~r/2, S2 )\T{J^ M (z/2)J^(-z/2)}\N(p + r/2, Sl )) . (12) 

The latter appears in the invariant matrix element (i.e. in the T-matrix) of the specific wDVCS process, and will 
be computed at the twist-2 level in the DVCS kinematics defined above. This approximation amounts to neglecting 
contributions of the order M 2 /q 2 and t/q 2 . Since the final state photon is on-shell, q 2 — 0, it follows in this particular 
kinematics that r • q± ~ —qf/2 — xb (pi ■ qi). Hence the momentum transfer should have a large component in the 
direction of the average nucleon momentum, 

r = 2r/p + A , (13) 

characterized by the skewness parameter 

_ r ■ q 



2 (p • q) 



(14) 



The remainder A in Eq. lfl3]] is transverse to both p and q [2l|, [22|. Moreover, in the DVCS kinematics, r) coincides 
with the scaling variable £ [2], where 



2 (p • q) 

It is then easy to verify that in the limit t/q 2 — > 0, one has £ = Xb/ (2 — Xb). 



q 2 

^ = ■ (15) 
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Having introduced the reduced wDVCS amplitude and the DVCS kinematics, we now turn to the formal light- 
cone expansion of the time-ordered product T {J^^iz/2) {— z/2)} in the coordinate representation. The expansion 
is performed in terms of QCD string operators [20J]. The string operators have gauge links along the straight line 
between the fields, however, for brevity we will not write them explicitly. The leading light-cone singularity is given by 
the s- and it-channel handbag diagrams shown in Fig. [H The hard part of each of the diagrams begins at zeroth order 
in a s with the purely tree level diagrams, in which the virtual weak boson and real photon interact with the (massless) 
quarks. The free quark propagator S(z) between the initial and final quark fields in the coordinate representation is 
given by: 



W - n : ^ - 7^4 e-*'«rhn ■ ( 16 ) 



2tt 2 (z 2 - i0) 2 J (2tt) 4 6 I 2 + iO 

Since the weak current couples to the quark fields through two types of vertices, qqZ° and qqW^, the quark fields 
at coordinates ±z/2 can carry either the same or different flavor quantum numbers. We shall treat these two cases 
separately. 






Figure 2: Handbag diagrams for the weak deeply virtual Compton scattering amplitude in the s-channel (left diagram) and 
u-channel (right diagram). 



B. Weak neutral amplitude 



In this section we expand the time-ordered product of the weak neutral and electromagnetic currents. Omitting the 
overall vertex factor — |e| g/ cosOy/, where g is the weak coupling constant, and 9w is the Weinberg angle, one has: 



iT 



K M (z/2) (-2/2)} = iJ2Q f $/(*/2)7"i P(*)-fl(4 -754)V/ [-z/2) 

+$f (-z/2) 7 "± (4 ~ 75^) % {-z) r^Pf {z/2) , 



(17) 



where Qf denotes the electric charge of the quark with flavor / (in units of |e|). In the Standard Model, the weak 
vector and axial vector charges are given by: 



u,c,t 



u y = 1/2 - 2<2 u , Cjt sin 9w and Cy S ' = -1/2 - 2Qd. s ,b sin 6\y > 
cY' 1 = 1/2 and c d /- b = -1/2 , 

with sin 2 dw — 0.23. Using the 7-matrix formula: 



(18) 
(19) 



(20) 



where s^ pvri = g pp g UTI + g m g pv — g^ 1, g pn is the symmetric tensor and e^ 1 " 7 the antisymmetric tensor [36j in Lorentz 
indices {n,v), the original bilocal quark operators can be expressed in terms of the vector and axial vector string 
operators with only one uncontracted Lorentz index: 



Of±(z\Q) = [i> f (z/2) Tq i; f (-z/2)±(z^-z)} , 
Ofi(z\0) ee [i, } {z/2) W i, f {-z/2)±{z^-z)] 



(21) 
(22) 
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Accordingly, the time-ordered product of currents in Eq. (T7j) assumes the form: 

iT{J£ M (z/2)Jfr N (-z/2)} = [s^Of-izW-ie^Qf+izlQ) 

77 Z f 

- c f A [s^O f ~ (z |0) - ie^Of+ (z |0)] } . (23) 



In contrast to the standard electromagnetic DVCS process, here we have two additional terms. Namely, the presence 
of the axial part 750^ of the V-A interaction gives rise to a vector current symmetric in {{i, v), and to an axial vector 
current antisymmetric in {fj,, v). 

The string operators in Eqs. ([2~T]) and (|22|) do not have a definite twist. To isolate their twist-2 parts, one uses a Taylor 
series_ expansion in the relative coordinate z. The expansion gives rise to the local operators ipf (0) 7r;-D Ml • ■ • D^ipf (0) 
and ipf (0) 7r;75-D M1 • • ■ D^ipf (0), where is the covariant derivative. To obtain the twist-2 contributions, one needs 
to project the totally symmetric, traceless parts of the coefficients in the expansion, which is effected by the following 
operation: 

[°f (* 1° )] syrn = ^ J & fa Mf ^^1?) T (z ^ -z)] , (24) 

°lt ( z 1° )1 = 7T1 I W fa ^ z / 2 ) (-^W T(z^-z)} . (25) 

J sym OZ I Jq 

The subtraction of traces is implemented by imposing the harmonic condition on the string operators on the right- 
hand-sides of Eqs. l(2~4|) and (|25l> . In other words, the contracted operators 

(z |0) = [$j (z/2) ti> f (-z/2) ± (z ^ -z)] , (26) 
0{± (*|0) = [j>, {z/2) M t {-z/2) ± (z ^ -z)] (27) 

should satisfy the d'Alembert equation with respect to z: 

d\ [0/±(*|O )] twiBt _ 2 = 0, (28) 

and similarly for the twist-2 part of G§ {z |0). 

To compute the amplitude in Eq. (fl2|) . we sandwich the contracted twist-2 operators between the initial and 
final nucleon states. To construct a parametrization for the nonforward nucleon matrix elements, we use a spectral 
representation, with the relevant spectral functions corresponding to off-forward parton distributions (OFPDs) [37I . 
Since the coordinate z runs over the whole four-dimensional space, in principle the parametrization should be valid 
everywhere in z. However, the inclusion of the z 2 terms in the matrix elements generates M 2 /q 2 and t/a 2 corrections 
to the amplitude (analogous to the well-known target mass corrections in deep inelastic scattering I23I 24l l25j|). and 
will hence be neglected. It is sufficient therefore to provide a parametrization only on the light-cone [22j, namely: 

<iV (P2, *a)l C» /=fc |0) |iV (pi, ai)>, a=0 = u(p 2 ,s 2 ) J dxe ixpz Hf{x,^t) 

+u{p 2 ,s 2 ) ^tJJl u^s^ J^dxe^ z Ef{x,U) , (29) 

(AT (p 2 , 52)10^(^0)1^^,5!)) „ = u{p 2 , S2 ) ^ 5 u{ Pl , Sl ) [ dxe^ z Hf{x^,t) 



-u{p 2 ,s 2 ) { -^ l5 u{ Pl , Sl ) J^dxe^ z Ej{x,U) ■ (30) 

The flavor dependent OFPDs in Eqs. ([29)) and lf30j) refer to the corresponding quark flavor / in the nucleon N. They 
depend on the usual light-cone momentum fraction x, the skewness parameter £, which specifies the longitudinal 
momentum asymmetry, and the invariant momentum transfer t to the target. As illustrated, for example in the s- 
channel diagram of Fig. [2j the parton taken out of the parent nucleon at the space-time point "-z/2" carries a fraction 
x + £ of the average nucleon momentum p, while the momentum of the reabsorbed parton at the space-time point 
"z/2" is (x- OP- 
Note that 05 ± {z\0) has a superscript opposite in sign with respect to the corresponding "tilded" distributions 
Hf ( x i€it) an d Ej {x,£,t). While the standard DVCS process gives access only to the plus distributions (i.e. the 
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sum of quark and antiquark distributions), scattering via the virtual weak boson exchange probes also the minus 
distributions. The latter correspond to the difference in quark and antiquark distributions, or the valence configuration. 
The advantage of using the plus and minus distributions, as opposed to the usual OFPDs Hf (x,£,t), Ef (x,£,t), 
Hf (x,£,t) and Ef (x,£,t), which parametrize the matrix elements of operators ipf (z/2) /:ipf (— z/2) and ipf (z/2) / 
zjsipf (—z/2), is that they both have well-defined symmetry properties with respect to the scaling variable x. By 
transforming z — > — z and x — > —x, one can readily establish the following crossing symmetry relations: 



Hf(x) = THf(-x). 



Ef (x) 



TEf (-x) , 



Hf(x) = ±Hf(-,- i 
Ef (x) 



f 

±Ef {-x) . 



(31) 
(32) 
(33) 
(34) 



Substituting the parametrizations (|29)) and (|30|) into the right-hand-sides of Eqs. (|24|) and (|25|) . one isolates the 
twist-2 terms by taking the derivative with respect to z, and integrating by parts over the parameter /?, keeping only 
the surface terms with the arguments ipf (±z/2) and -0/ (± z /2)- Finally, the integral over z is carried out with the 
help of the inversion formula for $ (z): 



d"z e a - 



ihz 



L 



2tt 2 (z 2 - iO) 2 (I 2 + iO) 



(35) 



Here the momentum I is given by I = (xp + q), so that I 2 in the denominator of Eq. (|35"1) becomes I 2 — 2 (p ■ q) (x — £). 
The expression for the reduced weak neutral virtual Compton scattering amplitude in the leading-twist approximation 
can then be written (we will implicitly deal with twist-2 amplitudes henceforth): 



1 WN — 



4(p 



dx 



-ie 
.J 



f 

^P vr >l B 



_ x (x-Z + lO) 

u{p2,s 2 )l v u(p 1 ,s 1 )Hf (x,£,t) + u(p 2 ,s 2 ) v v u(p 1 ,s 1 )Ef (x,£,t) 



+ie fJ - upr >l p 



u{p2,s 2 )l 11 l 5 u{pi,Si)Hf {x,£,t) -u(p 2 ,s 2 ) ^-i b u{pi,si)Ef {x,i,i) 
(P2,s 2 )'7r 1 75u(pi,s 1 )Hj (x,£,t) -u(p 2 ,s 2 ) 7^75« (Pi, «l) EJ {x,£,t) 

u(p 2 ,s 2 )i v u(pi,sx)HJ {x,£,t) + u(p 2 ,s 2 ) ^ L —T—-^-u(pi,s 1 )Ej (x,£,t) 



(36) 



This can be further simplified by using the light-cone (Sudakov) decomposition of the 7-matrix: 



7^ = fix + V ii 2 + 7± , 



(37) 



where the four-vectors m and n>2 in Eq. I|37p are light-like, n 2 



n 2 = 0, and satisfy the condition n± ■ n 2 = 1. 
Identifying n\ <-> p and n 2 «-» q 2 , and neglecting the transverse component 7^ (since it corresponds to the higher-twist 
contributions), Eq. J37|) takes the form: 



r 



(P ■ Q2, 



(?2 M 



p« 42) 



(38) 



Now, using the above decomposition, the Dirac equation, tfiu(pi,s\) =^2^(^2,52) = (recall that we neglect the 
nucleon mass), and the symmetry properties in Eqs. f31j) - lj34|) . the amplitude (|36|) becomes: 



1 wn — 



1 



4 (p • q) 



(p ■ Q2) 



(£i *) u (P2i s 2) 42U (pi, Si) + £w N (£, t) U (p 2 ,S 2 ) ■ 



42 j- j 42) 
AM 



u (pi,s a 



— _ — _ (02 ' V) _ 

-KwN (£,t)u(p2:S 2 ) 4275U(pi,Sl)+8w N (£,*) 2M U (P2,S2)l5U(Pl,Si) 
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where 



(p ■ 92) 



n WN (£> l ) U (P2) s 2) ^275"" (Pi, Si) - £^ N (£, t) 



(g2 • r) 
2M 



~ n WN (£> *) U (P2i s 2) ^2" bl, Si) - f^jv (£, *) U (p 2 , S 2 ) 



4M 



u(pi,Si) 



dx 



H 



_i(s-£ + iO) f 



+(-) 



(a, CO 



(x-£ + iO) / 

g?x i?/ (a;, £, t) 

da; ~ 
.! (z-f + iO) J 

e?x 77/ (a;, £, t) 

fia; 



£+< >(a;,£,f) 



j« «/ — 1 

I^ ( 3/4(A)/ i 



± 



1 



£ + iO x + £ - iO 



± 



1 



x - £ + iO x + £ - iO 

1 1 

x - £ + iO T x + £ - iO 



( -f+(-)(i,C,t) 

[ac — S + tOJ 7 



g?x -E/ (x, £, t) 



1 



-7T T 



x-£ + iO x + £-iO 



(39) 



(40) 



(41) 



(42) 



(43) 



is a new set of functions given by the integrals of OFPDs. The amplitude T^ N has both real and imaginary parts, 
with the real part obtained using the principal value prescription, and the imaginary part coming from the singularity 
of the expression 1/ (x =p £ ± iO). The latter generates the function 5 (x =p £), which constrains evaluating OFPDs at 
the specific point x = ±£. 



C. Weak charged amplitude 

For the weak charged amplitude, we write the expansion of the time-ordered product of the weak charged and 
electromagnetic currents (again omitting the overall vertex factor — |e| g/Vty as: 

iT{J EM {z/2)r wc (-z/2)} = --^ ^[Q f ^( Z /2) T Yf(l- 75) *f(-z/2) 

-Q / ^(-z/2) 7 "(l- 75 ) 7 " 7 ^/b/2)] • (44) 
Using the 7 -matrix decomposition in Eq. (|20|) . and collecting terms, one has: [H] 

iT{r EM (z/2) r wc (-z/2)} = -zfeE { s ^ n [Qf'$r (*/ 2 )^/ (-V2) -<?/(*-> -*)] 

/,/' 



-ie"""> [Q/^VV («/2) 7,75^/ (-«/2) + Q/ -» -*)] 
-s^" [Q,^ (z/2) 7 „ 75 V>/ (-z/2) -Q f (z^> -z)] 

■fie"'"" [Q f 4f> (z/2) 7 „V/ (-z/2) + Q, (z - -z)] } . (45) 
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Here the sum over quark flavors is subject to an additional condition, Qf — Qf — ±1, due to the fact that the virtual 
boson W carries an electric charge, so that the initial and final nucleons correspond to different particles. In other 
words, the neutron-to-proton transition is via the exchange of a W + using a neutrino beam, while the proton-to- 
neutron transition is via W~ using antineutrinos. Moreover, the vector and axial vector string operators are not 
diagonal in quark flavor, and are also accompanied by different electric charges. We can express the corresponding 
contracted string operators, which appear when extracting the twist-2 part of the expansion |45|) . as the linear 
combinations: 



of the operators: 



[Qf4f 0/2) jhpf (-z/2) ±Q f (z-> -z)] = Q ± Of'f+{z\Q) + Q T Of'f-{z\Q) , 
[Q f 4f, (z/2) M f (-z/2) ± Q/ (z - -z)] = Q ± O f 5 ' f+ (z\0) + Q T o( f - (z |0) 

0^±(z|O) ee [$ f , (z/2) ;% (-z/2) ± (z- -z)\ , 



O 



f'f± 



(z|0) = $ f ,(z/2) M f {-z/2)±{z^-z)] 



(46) 
(47) 

(48) 
(49) 



with the coefficients Q± = (Q f ± Q /) /2. The matrix elements of these operators are then parametrized in terms of 
the flavor nondiagonal OFPDs: 



<iV'(p 2 , S2 )|O^ ± (z|0)|iV(p 1 , Sl )) z2=0 = u(p 2 ,s 2 ) / dxe^ z Hf lf {x,i,t) 



+u(p 2 ,S 2 ) 



t). 



4M 



(N' (pa, s 2 )\ O f 5 ' f± (z |0) \N ( Pl , Sl )) z2=0 = 



: ,(p 2 ,s 2 ) fr B u(pi,si) J dx e lxp ' z Hj, f (x,£,t) 
-u( P2 ,s 2 )^- 75 u( Pl , Sl ) ( dx e^ z Ej, f (x,U) 



(50) 



(51) 



These correspond to the situation in which a quark with flavor / is removed from the target nucleon, and a quark 
with flavor /' is reabsorbed. Finally, with the help of Eqs. (|46|) and f47j) and the parametrizations (|50|) and (|5"Tj) . one 
obtains the twist-2 result for the reduced weak charged virtual Compton scattering amplitude: 



1 wc — 



1 



4 (p • q) 



1 



(P-Q2. 



u wc (£> *) u (P2, s 2 ) 4 2 u (pi, si) + £^ c (£, t) u (p 2 , s 2 ) 



4M 



u{p\,s x ) 



- n wc (Z> l ) u (Pa> s 2) ^275"" bi, si) + £ wc (£, i) %77^u (pa, s 2 ) 7 5 U (pi, Si) 



2M 



-ie^ pr > 



Q2pP v 



ft\VC (f » *) U (P2, S 2 ) ^275"" (PX, Si) - (£1 



(g2 • r) 
2M 



«(P2,S 2 )75M(Px,Sx) 



"W ffc (£,i)u(p 2 ,S 2 ) gf 2 w(px,Sx) (^ t ) M (P2,S 2 ) 

2 



4M 



-tt(pi,Si) 



(p- 92) 

Tx (t) u (p 2 , s 2 ) ^ 2 u (pi, si) + ^ 2 (t) u (p 2 , s 2 ) ^ 2 4^/ ~ M (Pi> s 
-0^ (i)w(p 2 ,s 2 ) ^ 2 7 5 u(pi,si) + £ P (i) g ^ (p 2 , 82) 75^ (px, sx) 



2M 



where the integrals of OFPDs are now given by: 

" x dx 



H 



/,/' 



x (1 - £ + »0) 



Q +( -)H+ f (x, £, t) + Q-{+)Hj, f (x, £, t) 



(52) 
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and 



E 

/,/' 

E 



Qf 



± 



Qf 



dx Hf t (x. £. t) , 
_ x ' M ; \x-£ + i0 x + £-i0 



dx 



Ax-£ + iQ) 



I dxE rj {x,t,t) 
fJ' J - 1 



Q+(-)E+ f {x, t, t) + Q- W Ej lf (x, £, t) 

Qf ± Qf 



/,/' 



E 



dx 



! (I - £ + *0) 



x-£ + z0 x + ^-iO 
Q+(-)Hp f (x, 1 1) + Q- (+) HJ, f (x, & i) 



! ; ; V ' S ' ; \x-£ + iO ' x + £-iO 



SwPiU) ^ E 



dx 



E 



Q+ ( -)E+ f {x,i,t) + Q- [+) Ej lf (x,i,t) 

Qr ^ Qf 



dx E f ,f{x,Z,t) I ^— — T , 

i \x — t; + iv x + £ — IV 



5p (t) 



Qf) / dxE flf (x,Z,t) 



J2Q- f dxE rf (x,tt) = J2(Qf-Qf) I 

E<2- / da; ^/'/ (*> ^ f ) = E (Qf - <w / da; Hft *) 
Eo- / dx Ef>f {x, e, = E (^/' - Q/) / dx ( x ' *) 



(53) 



(54) 



(55) 



(56) 

(57) 
(58) 
(59) 
(60) 



It is important to note the difference in the structure between the weak neutral and weak charged amplitudes. The 
latter has an additional term, which is symmetric in indices /z and v and is determined by the form factors lf57j) - 
(|60"ll . This completes the derivation of the wDVCS amplitudes. 



III. WEAK AND ELECTROMAGNETIC DVCS PROCESSES 

In this section, we consider specific examples of wDVCS processes for neutral and charged currents on an unpolarized 
nucleon target. For the former, we examine both neutrino and electron (or charged lepton in general) scattering, while 
the latter is illustrated using the neutron-to-proton transition in neutrino scattering. Before we examine the specific 
processes, however, we discuss the kinematics which are common to all DVCS-like reactions. 



A. Kinematics 



The generalized DVCS process: 

l(k)+N( Pl ) — > l'(k')+N'(p 2 )+^(q 2 ) , (61) 

where a lepton I (with four-momentum k) scatters from a nucleon N (pi) to a final state V (fc'), nucleon N' (j>2) and 
a real photon 7 ((72), is illustrated in Fig. [3j The wDVCS diagram, or Compton contribution, is depicted in Fig. Ufa), 
which corresponds to the emission of a real photon from the nucleon "blob". The other two diagrams, (b) and (c) 
in Fig. [3j illustrate the Bethe-Heitler process, where the real photon is emitted from either the initial or final lepton 
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leg. Here the nucleon "blob" represents the electroweak form factor, while the upper part of each diagram can be 
calculated exactly in QED. 




a be 



Figure 3: Weak DVCS (a) and Bethe-Heitler (b and c) diagrams contributing to the leptoproduction of a real photon. 

We denote the various four-momenta in the target rest frame as k — (lj, kj, p\ — ( M, qY k' = (uj 1 , k'J , P2 — 

(E 2 , P 2) and q 2 = (^2, 92)- The differential cross section for lepton scattering from a nucleon to a final state with a 
lepton, nucleon and a real photon is given by: 

a 1 irria 1 ru 1 u d 3 k'd 3 p 2 d 3 q 2 

where T represents the invariant matrix element containing both the Compton and Bethe-Heitler contributions: 

T = Tc + Tbh ■ (63) 
In the target rest frame, the total lepton-nucleon center of mass energy squared is s = (pi + k) 2 = 2Mw + M 2 , where 



we neglect lepton masses. Integrating the cross section in Eq. ([62|) over the photon momentum yields: 



da = ITI 2 ' 



4Mw (2tt)- 



(k+ Pl -k-p 2 )^ — . (64) 



i2 



The 5-function here provides the constraint s + M 2 — 2 [{v\ + M) E 2 — qi ■ P2] = 0, where the invariant s = (pi + qi 
and qi = k — k' = (i/i,qi) is the four-momentum of the virtual weak boson with the energy v\ — u> — u>' and the 
magnitude of the three-momentum \qi\ = \J v\ — q\ ■ In addition, we choose a coordinate system (depicted in Fig. [4j) 
where the virtual weak boson four-momentum has no transverse components, q\ = (^i,0, 0, \q\\), and the incoming 
and outgoing lepton four-momenta are k — u) (1 , sin <j>, 0, cos <f>) and k' — u>' (1, sin <fi' , 0, cos </>'), respectively. In this 
reference frame, the azimuthal angle of the recoil nucleon corresponds to the angle ip between the lepton and nucleon 
scattering planes. Using the (5-function in Eq. (|64|) to integrate over the polar angle of the outgoing nucleon, and 
expressing its energy in terms of the invariant momentum transfer t = (pi — p 2 ) 2 as E 2 = M — t/ (2M), we find: 

[\ + {u-u')IM\t + q\ 
cos 02 = 1^^^=^^^= , (65) 



2\p 2 \ J{u-L0'y~q 2 



where \p 2 \ = y/-t[l - tj (4M 2 )], and the cross section takes the form: 



1 — ,2 1 



da = ^-l T l /0 s4 / n =dw' dE 2 d (cos <t>')d<p . (66) 

32MW ( 2 ^) yfc co') 2 q\ 



Here the angle 4>' denotes the polar angle of the scattered lepton. 
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y 




Figure 4: Kinematics of the generalized DVCS process in the target rest frame. 

Instead of the kinematical variables (a/, E 2 , cos </>'), it is more convenient to express the differential cross section 
in terms of the variables {Qf,t, Xb), where Q\ = —qf. In the target rest frame, we have Q\ = 2Mu>yx B , where the 
invariant y = p\ ■ q\/pi ■ k = (uj — ui') /u> and accordingly, Eq. f65|) turns into: 

x B M (1 - t/Ql) - t_[ (2M) 

cos 02 = — ; — • 67 

\p 2 \ y/l + 4a? B M*/Qi 

Furthermore, from the conservation of three-momentum, q\ = k — k' , and writing the invariant Q\ = — (k — k') 2 = 
2ujlo' (1 — cos 4>i ep ), where <fii ep is the angle between the incident and scattered lepton momenta, we get for the polar 
angles of the incoming and scattered leptons: 

1 + 2yx B M 2 /Ql 



and 



1 



t 2yx%M 2 

{i-y)Q\ 



y/l + 4x* B MyQl 
respectively. The differential cross section can then be written as 

d 4 a 1 1 x B y 2 



(69) 



T| 2 . (70) 



dx B dQ(dtdip 32(2tt) 4 Qi ^l + ix 2 B M 2 /Ql 

Finally, energy conservation, M -\-v\ = E 2 + 1*2, implies that the energy of the outgoing real photon in the target rest 
frame is v 2 — (Q 2 + x B t) / (2Mx B ). Alternatively, we can write the invariant momentum transfer as t = (q-2 — qi) 2 — 

—Q\ — 2i>2 (uy — \Ju 2 y 2 + Q\ cosOb^J , where 6*b 7 denotes the scattering angle between the incoming virtual weak 
boson and outgoing real photon. Combining both expressions gives: 

Q 2 u 2 + x B M (Ql + t) 

COS0 B ~, = .. i (71 

or, on the other hand, the invariant t can be expressed as a function of the angle #b 7 : 



1 + 1 - v/i + 4x|M7Q2 cos Q Bl Q\l (2x 2 B M 2 ) 
1 1 + (l - y/1 + Ax B M 2 /Q\ coaOBy) Q\j (2x B M 2 ) 
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In the following, we determine the kinematically allowed region for the generalized DVCS process, which requires 
finding the upper and lower limits of the invariants xb, y and t. The kinematics are subject to the following constraints: 

1. The energy of the incoming lepton beam is fixed at u> = 20 GeV. 

2. The invariant mass squared of the virtual boson-nucleon system should be above the resonance region, s = 
(pi+<Zi) 2 >.Smi„ = 4GeV 2 . 

3. The virtuality of the incoming boson has to be large enough to ensure light-cone dominance of the scattering 
process, Q\ > Q 2 lmin = 2.5 GeV 2 . 

4. The momentum transfer squared should be as small as possible, e.g. 0.1 GeV 2 < \t\ < 0.2 GeV 2 , which yields a 
low-energy nucleon and a high-energy real photon with in the final state. 

The constraint (2) leads to a lower limit on the variable y: 

_ s mia - M 2 

Vmin ~ (s — M 2 ) {1 — xb) ' [ ' 

On the other hand, when the incoming lepton is aligned along the z-axis, i.e. <pi ep = 180°, the scaling variable y 
reaches its maximum value: 

ymax =(l + ^J ■ (74) 

The region in the xsy plane is then bounded by three curves, illustrated in Fig.[5j These are given by y m ; n in Eq. (|73| . 
Z/max in Eq. (f74|) . and y = Q 2 min / [(s — M 2 ) xb], which follows from the constraint (3). Next, both the lower and 
upper limits of the invariant t can be determined in the virtual weak boson-nucleon center of mass frame. Here the 
kinematic limits of t are given by (up to relative corrections of the order xbM 2 /Q\): 



-M 2 x 2 B _ M 2 x 2 B - 2M 2 x B + Q\ {x B - 1) (xb 

tmin 1 xb (1 M 2 IQX) and tmax _ 1 - x B (1 - M 2 /Q\) ' (?5j 



at the scattering angles of 0° and 180°, respectively, between the initial and final nucleons. 

y 



0.8 



0.6 - 



/ 

/ 

/ 



/ 

/ 

/ 



0.4- . / 



0.2 0.4 0.6 0.8 



Figure 5: Kinematics of the DVCS process: the dashed curve gives j/ m i n in Eq. I|73p . the solid represents y max in Eq. I|74p . and 
the dotted corresponds to y = Qi m i n / [(s — M 2 ) xb\ ■ The region bounded by the three curves is that kinematically allowed 



for s > 4 GeV 2 and Q\ > 2.5 GeV 2 with an lo = 20 GeV lepton beam. 



Finally, we select one kinematical point within the allowed region in the xsy plane, namely Q\ = 2.5 GeV 2 and 
xb = 0.35. For this point, in Fig. [6] we plot the invariant momentum transfer against the scattering angle 6*b 7 , which 
varies between —0.15 GeV and —1.433 GeV for < 0B-y < 15°. We also set the angle between the lepton and 
nucleon scattering planes to (p = 0. Then, with the so-called in-plane kinematics, the polar angles of both incoming 
and scattered leptons are fixed to 4> = 20.2° and (/)' — 25.3°, respectively. 



14 




Figure 6: Invariant momentum transfer squared as a function of the scattering angle between the incoming virtual weak boson 
and outgoing real photon in the target rest frame, for Q\ = 2.5 GeV 2 and xb = 0.35. 



B. Model 

To proceed with a numerical study of the wDVCS process, we need a model of the OFPDs. Since this is an 
exploratory study, rather than a detailed one intended for comparison to experiment, we choose for simplicity a toy 
model to illustrate the main features of wDVCS and the differences with electromagnetic DVCS. 

Inherent in our model are the following approximations: 

• The sea quark contributions are negligible, which implies that the "plus" OFPDs coincide with the "minus" 
OFPDs. For that reason, Hf = HJ al and Hf = HJ al for quark flavor / = u or d, and similarly for the Ef and 

Ef distributions. 

• The t dependence can be factorized from the other two scaling variables (x and £) for all distributions, and the 
dependence of the OFPDs on t is characterized by the corresponding form factors. 

• The £ dependence of OFPDs appears only in the Ef distribution. 

The parametrization of the unpolarized quark OFPDs in the proton is taken from Ref. [2(|: 

H™ l (x,£,t) = u; al (x)F lu (t)/2, (76) 

Hr l (x,U) = d; al (x)F ld (t) , (77) 

E%*(x,t,t) = u v p al (x)F 2u (t)/2, (78) 

E v d al (x,U) = d; al (x)F 2d (t) , (79) 
where the unpolarized valence quark distributions are given by [271 ]: 

u v p al (x) = l.89x-° A (1- a;) 3 ' 5 (1 + 6x) , (80) 

d v p al {x) = 0.54a;- - 6 (1 - xf 2 (1 + 8s) . (81) 

They closely reproduce the corresponding GRV parametrizations [28[ at a low normalization point, — qf ~ 1 GeV 2 
[111 . The u- and c?-quark form factors in Eqs. f76j) - ([79"ll can be extracted from the proton and neutron Dirac and 
Pauli form factors, neglecting strangeness and heavier flavors, according to F lp ^ = Q u F\ u (d) + Q dF\d( u ) and ^2p(n) = 
QuF 2 u(d) + Qd,F2d( u ) j which are in turn related to the Sachs electric and magnetic form factors: 



G E p(n) (t) - J^G M p(n) {t) 



F lp ( n ) (t) = 

F 2 p(n) (*) = [G M p(n) (t) ~ G E p(n) (*)] fl - ) ' ( 83 ) 
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At small t, both the proton and neutron magnetic form factors, as well as the proton electric form factor, are 
approximated by dipole fits, while the electric neutron form factor is taken to be zero: 



Ge p (t) 



Gm p (t) _ Gmu (t) 

1 ~t~ K n 



t 

A 2 



and G En (t) = , 



(84) 



where n p = 1.793 and K n = 
parameter A 2 = 0.71 GeV 2 . 



- 1 .913 are the proton and neutron anomalous magnetic moments, respectively, and the 



In the polarized case, we take for the valence distributions the factorized ansatz from Ref. [30|: 

-2 



TTl A 



t 



Hr i (x,a,t) = Ad; ai (x)[i- — 

F > mil 



(85) 
(86) 



where the mass itia = 1-03 GeV. The polarized valence quark distributions in the proton can be expressed in terms 
of the unpolarized distributions using a parameterization motivated by the SU(6) quark model [3lL l32|: 



A vat ri I val ^ ival 

U p = COSUd \ U p — -Cl p 

1 



(87) 



Ad v p aL = COS0D 



_jval 

3 p 



where cos6»d = [l + 0.06 (l - x 2 ) 
is dominated by the pion pole: 



Finally, for the Ef distribution we use an ansatz in which the distribution 



1 



'(N<0 



x + £ 



2 " w 2^ 

-K al (x,u) ■ 

The function (t) is assumed to have a form valid for —t <C M 2 [33]: 

1.7 (ml - t) /GeV 2 



(m 2 - t) 



1 - 



(1 - 1/2 GeV 2 ) 2 



(89) 
(90) 

(91) 



where m w ~ 0.14 GeV, and the axial charge of the nucleon is 17,4 ~ 1.267. For the pion distribution amplitude in 
Eq. |90|) we choose, for simplicity, its asymptotic form: 



(j)^ (it) = 6u (1 — u) 



(92) 



Having defined the kinematics and described the model for the nucleon OFPDs, in the rest of this section we 
examine several specific DVCS processes. 



C. Weak neutral current scattering 

In this subsection we discuss two examples of the weak neutral current scattering process, for neutrino-proton and 
electron-proton scattering. For the case of an unpolarized neutron target, one can use isospin symmetry to express the 
neutron OFPDs in terms of proton OFPDs. To make our comparisons meaningful, in all cases the relevant unpolarized 
differential cross sections f70|) will be plotted as a function of the angle 6*b 7 for the same kinematical point, namely, 
Q\ — 2.5 GeV 2 and xb — 0.35, with a lepton beam energy uj = 20 GeV. 



1. Neutrino-proton scattering 



Since neutrinos do not interact with photons, neutrino scattering from a proton via the exchange of a Z° boson 
measures the pure Compton contribution, with no contribution from the Bethe-Heitler process. The T-matrix for this 
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process 
iT 



Up 



u(k') 



-if] 



cos 6w 



7 



u(Jfe) 



j (.gi^A - glt/glA /M|o) 



COS 6*iy 



< fe) , (93) 



is then solely given by the DVCS diagram, Fig. OJa). The vector and axial vector couplings at the vvZ° vertex 
are c v = c v A = 1/2. By taking into account that M§ S> Q 2 and cosOw = Mw/M z o, and further recalling that 
g 2 / {d>M w ) = Gf/V%, where Gf — 1.166 • 10~ 5 GeV -2 is the Fermi coupling constant, we arrive at: 

T„ p = V2 \e\ G F u (k') 7 , (1 - 75) u(fc)e* (92) . (94) 

To compute the unpolarized cross section, one should average the square of the amplitude T„ p over the spins of the 
initial particles, and further sum it over the spins and polarizations of the final particles. In particular, the summation 
over the photon polarizations is performed using the Feynman gauge prescription, i.e. one can replace 



(82) Cc. (©) 

7 polar. 



(95) 



by virtue of the Ward identity. As a result, we then write the spin- averaged square of the T-matrix in terms of 
factorized neutrino and weak neutral hadronic tensors, 



IT^I 2 = 8naG 2 F L%Htf N , 



(96) 



(97) 



where a = e 2 / (4tt) ~ 1/137 is the electromagnetic fine structure constant. The neutrino tensor L„p is given by 

L ip = 8 [kvk'fj + kpk' v - g v p (k ■ k') - ie v ^ T k a k' r \ . 

The weak neutral hadronic tensor H^ N , on the other hand, has a considerably more complicated structure. In the 
DVCS kinematics, however, in which 0{t/q 2 ) and (D(M 2 /q 2 ) terms are neglected, and using the Feynman gauge 
prescription, this reduces to: 



H 



up 
WN 



1 



-p 



2 *WN y^iiWN 

1 



WN 



up 



u P 

Q212 



■-2WN 



(P- 92) 



ie^ 5X p 5 q 2 x 



(98) 



where the functions C\wn and C^wn are given by: 



WN 



(i-e) [\n 



i+ 1 

l WN\ 



n 



WN 



H 



WN 



n 



WN 



-u 



AM 2 ) 



'WN I 



'WN I 



AM 2 



°WN 



'WN 



2£ 3? (h wn £ jy N + T~l WN £ WN + T^-wn^wn + ^wn^wn) > 



C 



2WN 



n WN n WN n WN n WN 



AM 2 



^ (Z-wn^wn) € 



AM 2 



( 



3^ I £w/ 



£ ^ y^WN^WN + ^WN^WN + ^WJV^WJV + ^WN^WN J 

The result of the contraction of the tensor H^ N with L^l in Eq. ([97"1) can be presented in a compact form as: 
_ \<o-KaG F Q 2 



T 



Ci 



WAT 



l-(l-tf)' 



c 



2WN 



(100) 



(101) 
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2. Electron-proton scattering 



Unlike in the neutrino case, for electron scattering both the Compton and Bethe-Heitler processes contribute, in 
which case the ep amplitude is given by the sum T ep = Tc ep + T^BHep- By replacing neutrinos with electrons in 
Eq. lj93j) , we can immediately write for the Compton part of the T-matrix: 



T Cep = 2V2 \e\ G F u (A') 7, (c e v - 75 ci) u(A)e* (<fe) , (102) 
where the couplings are now Cy = — 1/2 + 2sin 2 9w and c\ — —1/2. The spin-averaged square of Tc e p then reads: 

(103) 



|T Cep | 2 = 32iraG 2 F L [ ^H^ N 



The hadronic tensor H^ N is the same as in Eq. J98|), while the electron tensor has an additional factor of 1/2 from 
averaging over the initial electron spins: 



r(e) _ 9 



{ ( c y) 2 + ( c a) 2 [kuk' f j + kpkl - g„f) (k ■ k')] ~ 2c v c e A ie u fjark (J k' T ^ , 



(104) 



and hence 



-Cep\ 



^ _ 16iraG F Qj 
V 2 



{ [{cvf + (ca) 2 ] [l + (1 - vf] Ciwn - 2c e v c A [l - (1 - yf\ C 2W n) ■ (105) 



For the Bethe-Heitler contribution, since both the initial and final leptons are electrons, both diagrams (b) and (c) 
in Fig. [3] contribute. The full Bethe-Heitler amplitude is given by: 



7 M (l'+ 42) I" (c c y-l5C A ) , 7"(<fr-75<£)(ff- h)^ 



(k' + q2Y 



T BHep = 2V2\e\G F e;{q 2 )u(k') 

X (P(P2,S 2 )\J? C (0)\ P ( P1 , S1 )) . 

Accordingly, the electron and hadronic tensors in the Bethe-Heitler amplitude 



u(k) 



are given by: 



-W3 



1 



|Ts_ffe P | 2 — 32iraG F L U g H H^ 



(7" ^2 + 2fc'") 7" (c e y - 75ci) , 7" K - 754) (^7 M " 2fc") 



(106) 



(107) 



2(fc'-g 2 ) 2(fc.«ft) 
7^ - 75^) + 2*;) , ( 7m ^2 - 2k,) ^ (cf v - lb c A ) 



(A' • q 2 ) 
2 



+ 



2 (A' • ffl) 2 (A • q 2 ) 

{ [(cVf + (c A ) 2 } \W<& + k'W - g vti (A ■ q 2 )] - 2c^ie^ ffT A CT(?2T } 
{ [(c e vf + [c A f] [k^ql + ~ 9 V0 (k' ■ q 2 )] + 2c e v c e A^ Uf3aT Kq2r} 



(A • q 2 ) 

— 1- { \(c c v f + (c A f] [(A ■ q 2 ) [Wk 113 + k 13 ^ + 2 W] - (k' ■ q 2 ) [k v k l!i + k^k lv + 2A"A /3 1 

(A' ■ q 2 ) (A ■ q 2 ) L L J 

+ (k ■ k') \k v 4 + k P q 2 - k w ql - k'^q" + 2k v k' (i + 2A /3 A" y l + 2g vf} [k ■ k') [(*' ■ g 2 ) - (A • g 2 ) - (A ■ k')] 



and 



-2^^^^ [(A • A') [2k a k' T + k a q 2T + k' a q 2T \ + [(A ■ q 2 ) - (A' • «,)] A CT A' T ] | 



H ™ = 5E<p(P2. s 3)l^ c ( )b(Pl.*l))<p(P2,^)l^ c (o)|p(Pl > *l)) , 



(108) 



(109) 
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respectively. The weak neutral transition current J^ c in Eq. (|109p consists of vector and axial vector parts: 

J? c (0) = \ [K NC (0) - < c (0)] 

= ^J2[ c v^ ( )7^/(0) -c^/(0) 7,75^/(0)" . 



(110) 



Their matrix elements are parametrized as (see e.g. Ref. [34]): 



{p( P2 ,s 2 )\V^ c (0)\p( Pl , Sl )) = u(p 2 ,s 2 ) 



JVC fj.\ 1Vu)T 



2M 



u(pi,s 1 



(p (pa, S2 )| ^ c (0) |p (p x , Sl )) = u (pa, S2 ) [G^ c (t) 7,75 - G£ c (t) ^] u (p x , 5l ) 



(111) 
(112) 



which then gives the following expression for the hadronic tensor: 



H, 



BH 

V0 



9vP 



[F? c {t) + F» c {t)) 2 + 1 



(Ff c (t)Y 



AM 2 



(F 2 NC (t)Y 



+ 



(piu - y) (pifi 



2 



M 2 I 1 



4M 2 



9ufi 



NC u\\ 2 i r v r 



(g n a c (t)) 



AM 2 



■^p [GT (t) + G N P C (t)] 2 - ie^ffi K c (t) + F 2 " c (t)] GT it) 



(113) 



We express the result of the contraction of Eq. (| 1 13|) with the leptonic tensor in Eq. (| 108|) in terms of the kinematical 
invariants xb, y, Q\ and t, and the scalar product: 



k ■ q 2 — ujv 2 [I — (sui 4> sin 8 By + cos cf) cos 9 B-y 
1 + x B t/Ql 



Ql 



y(l + Ax 2 B M 2 /Q 2 ) 



i _ I _ I ( i + 2 ^l^ 2 

2 2 



l + t/gf 
l + a; B t/Q 2 



A 



i-y 



y 2 x 2 B M 2 



1 + t/Q? x|M 2 / 1 + t/Q 2 



l + x B t/Q 2 Q\ \l + x B t/Q 2 



(114) 



Moreover, it is convenient to introduce the dimensionless variables t = t/Q\, fi = M 2 jQ\ and k = (k ■ q 2 ) /Q\. The 
latter can be written in the following form: 



2y(l + Ax B p) {1 - 2/C - r [1 - , B (2 - y) + 2x%y,] + 2x%y,} - \ 



where 



K 



- (1 - y - x 2 B y 2 p) (1 - x B ) t + x b (1- x B ) t 2 + x 2 B y 2 ^ (1 + r)' 



(115) 



(116) 



is the l/yj Qf-power suppressed kinematical factor. In terms of these variables, the Bethe-Heitler squared amplitude 
can then be written: 



iTsjjepT = — 327raG 2 ? <2 2 



1 



2k (2k + t + l)x 2 B y 2 



(4) + (c e A y 



2 2 

x B y T 



(2k + l) 2 + (2k + t) 2 \ [F? c (t) + F 2 NC (t)Y 
+2 [2r(l-y) - x B yr (4k + r + 1) + y 2 [r + x s r (2k + r) + z|/z [l + 8k 2 + t 2 + 4k (1 + r)]]] 
x "« C W) 2 -4^p(^ C W) 2 ]" 

+ [4t - 2yr [2 + i B (4k + r + 1)] + y 2 [2r + 2a; B r (2k + r) + x% (t - 2p) [l + 8k 2 + t 2 + 4k (1 + r)]]] 

x "(4) 2 + (4) 2 ] (Gf W) 2 

+4.t b2/ t [2 (1 - r) + y [2 (2k + r) - x B (1 - r) (4k + r + 1)]] [ff c (t) + F? c (*)] G^ c (t) 1 , (117) 
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where the term 1/ [2k (2k + r + 1)] comes from the lepton propagators. Note that, in contrast to the Compton 
contribution, Eq. I|117p is the result of the exact calculation. Each of the vector form factors can be further decomposed 
into linear combinations of flavor triplet, octet and singlet form factors: 



Fi N (2)(t) = (1-2 sin 2 I 



(*) + g^i 8 (2) (*) 



6^°(2) (*) > 



(118) 



while the axial and pseudoscalar form factors can be written, in general, as differences between the isovector and 
strangeness form factors: 



1 



g a(p) (*) — G%(P) (*) ~ 2^( p ) ^ ' 
If one further neglects sea quark contributions, then these form factors can be written as: 



*1(3) (*) 



2 [ F l(2)u (*) - Fl(2)d (*)] 

*i(2) (*) = *i (a) (*) = [F H 2)u(t)+F md (t)] , 

G A(P) (*) = G \(P) (*) • 



(119) 

(120) 
(121) 
(122) 



Recall that the t dependence of the u- and d-quark form factors is given by the Dirac and Pauli form factors, whereas 
for the axial and pseudoscalar form factors we use the parametrizations: 



r<NC 



(t) 



9A 

2 



1 



t 



G% c (t) ( AM 



G» c {t) 



The interference terms between the Compton and Bethe-Heitler contributions, 

ep — -L Cep -L BHep + 1 Cep 1 BHep ; 



(123) 
(124) 

(125) 



are particularly interesting, since they are linear in the integrals of OFPDs. Substituting the expressions in Eqs. (|102p 
and (|106p into Eq. I|125p . and averaging and summing over the initial and final spins, respectively, the interference 
term can be written in terms of the electron and hadronic traces as: 



--ep 



-AnaG 2 F \ Tr \ }t' lv (c e v - l5 c e A ) ]6 
xTr 

+Tr I ^ 7(3 {c<y - l5 c e A ) ft' 



1 P (c e v ~ l 5 c e A ) + 2J£) ( 7|t fa - 2fc M ) ^ (4 - l5 c A ) 



2 {k> ■ q 2 ) 



+ M)T^ N (^ + M) 



f nc 



It7f3 7 



2M 



2(k-q 2 ) 

NC U s 75^/3 



2M 



( 7 " 42 + 2k'») r {c e v - 75^) , Y K - l5 c A ) (far - 2fc^) 



2 [k' ■ q 2 ) 



2 (k ■ q 2 ) 



xTr 



(A+M) (f[! WN )\fc + M) 



Fr(t) lv -F™(t). 



NC u\ 75>"i/ 



~2M A ' ' 1 15 p 



2M 



where the amplitude Ty^ N is the spinorless part of the reduced virtual Compton scattering amplitude 



WN- 



T wn = u(p2,s 2 )T^ N u(pi,si) 



(126) 



(127) 



To be consistent, we need to keep the same level of accuracy as in the Compton part, i.e. we should neglect terms that 
are of O (t/qf) and O (M 2 /g 2 ) order. Accordingly, the variable K becomes K = [(1 — y) — 2JC] / (2y). Furthermore, 
we should only keep terms linear in JC. After the contraction, we obtain: 



T 



1 



2k(2k + t + l)x B y 3 



{[(cv) 2 + (c A ) 2 } [4/C(2-2 y + y 2 )] 



F? c (t) n {H+ N ) - j^F 2 NC (t) n (S+ N ) + G N A C (t) n (H WN ) 



20 



[ Fl NC (t) + F 2 NC (*)] 5ft («+ + G N A C {t) [3? 



+44 [8/Cy (2 - y)] 



4A/ 2 



F 2 ivo (t)ft(£ 



G 



iVC 



(*)] B + (i) [H (W+ JV ) + K (£ + „)] 

D. Weak charged current scattering 



}■ 



(128) 



In the weak charged current sector, we consider neutrino scattering from a neutron via the exchange of a W + boson, 
producing a proton in the final state. The T-matrix of the Compton contribution in this case is: 



T Cl/n = V2 \e\ G F u (k') 7 , (1 - 7s) u (k) e* (q 2 ) T& v c , 



(129) 



where the amplitude T w u c ' ls > m g enera li given by Eq. lj52|) . In our simple model, however, the quark flavors / and 
/' in Eqs. (|53| - lj56|) are d and u, respectively, and the coefficients equal to Q+ — 1/6 and Q_ = 1/2. To proceed, 
we relate the flavor nondiagonal OFPDs to flavor diagonal ones using isospin symmetry. Specifically, the nucleon 
matrix elements (p{p 2 ,s 2 )\O ud± (z |0) \n (pi, si)) and (p (p 2 , s 2 )\ 0% d± (z\0)\n(pi,si)) in Eqs. J3S| and (|49j) that are 
nondiagonal in quark flavor are expressed in terms of the flavor diagonal matrix elements according to Ref. [35]: 

(p (pa, * 2 )| <T d± (z |0) |n (pi, si)) = (p (pa, s 2 )\ O u± (z |0) |p (p l5 Sl )) - (p (pa, a a )| O d± {z |0) |p (px, , (130) 

and similarly for (p (p 2 , s 2 )| Og d± (z |0) |n (pi, Si)). The convolution integrals in Eqs. (53]) - O and Eqs. |57j) - ([60]) 
then become: 



a H'C 


&*) = 


y da; 


(x,£,t) 


-H d (x,£,t)] 




(£,*) = 




(x,£,t)- 


-E d (x,£,t)] 




(£,*) = 






(x,£,t) 


-Hd(x,£,t) 


?+(-) 


(£,*) = 




E u 


(x,£,t) 


-Ed(x,Z,t) 



Qu 



± 



x-£,+i0 x + £,-i0 

Qu _j_ Qd 

x-£ + i0 x + £-i0 
Qd 



x- £ + i0 x + £-i0 

Qu Qd 

x-£ + iO T x + £-iQ 



and 



Ti (t) = (Qu - Qd) [Fiu (t) - Fi d («)] , 

T 2 (t) = (Q„ - Q d ) [F 2u (t) - F 2d (t)] , 

(t) = (Q u - Qd) <7A (*) , 

Gp (t) = (Q„ - Qd) .9p (t) . 

The spin-averaged square of the T-matrix in Eq. (| 129|) can then be written as: 



where the weak charged hadronic tensor is given by: 



jTTCXLrp h v p tl wc , 



(131) 
(132) 
(133) 
(134) 



(135) 
(136) 
(137) 
(138) 



(139) 



V0 

wc 



T^ v ( n~P \ * 
2 wc \ 2 fj,wc J 

\ { Ciwc 

2M 2 
(p-fe) 



1 - 



(P- 92) 



4M 2 



V 1 1/ fl 

p q 2 +p q 2 



C3WC p v q 2 + C 4 wc p p q 2 + 2C 5W c p v p 



M 2 



(P ■ «2) 



1 - 



AM 2 



V 

q 2 q 2 



-2WC 



ufiSX 



psq 2 \ 



(140) 
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with the following set of functions: 



Ci 



wc 



-e 



AM 2 



{\^-wc 


2 

+ 


I c wc 


2 

+ 



n 



- i 2 



n 



wc 



n 



wc 



! )-(« !+ 4i^)(K-| S + l^0 



'WC 



2£ 2 5? (U+ 



--2WC 



= -2 



--3WC = 



-AWC 



1 £ 2 ) y~L wc TL wc + T^wc^-wc) \£ 1 ^ 2 
— £ 5i (t^wc^wc + ^Trc^vrc + ^ifc^ivc + ^vrc^vrc) 

{i-e)(^np c + g A n w * c )-(e 



wc^vvc + ^ivc^'c + ^vrc^vrc + ^vrc^vrc) ' (1^1) 

\ /^2 * \ m t c + * c- \ /-2 t <rs ( o+* c- \ 



^ i^wc^wc) ^ 4^2 ^ (^VrC^VrC J 



4M 2 



2 C WC) 



~^lld 2 ~ ip pE wc) - ¥ (Fi^wc + ^ H wc + Sa£ w * c + GpHwc) > 



(1 - £ ) + 
t 



4M 2 



(•^2^0) 



5 1 re 



(1 - e) {Tl + Q\) U 2 + ^5) -ff - e 2 ^^ - 2^ 2 (V2 + SaGf) • 



(142) 

(143) 

(144) 
(145) 



As opposed to the weak neutral case, the additional functions C3WC, Ciwc and C$wc appear as a result of the extra 
term in the weak charged amplitude |52|) . Moreover, by neglecting the mass of the outgoing charged lepton, the 
leptonic tensor L^l coincides with the neutrino tensor L^l in Eq. l(97|). Consequently, after ignoring terms O (t/q 2 ) 
and O (M 2 /q 2 ), Eq. {139]) simplifies into: 



l6naG 2 F Q 2 



|Tc 



i { [l + (1 - y) 2 ] Cxwc - [l - (1 - yf] C 2WC } 



(146) 



For the Bethe-Heitler background, only diagram (b) of Fig. [3] contributes, for which the amplitude is given by: 



V2\e\G F e*(q 2 )u(k>) 



42) 7" (1-75) 



u(k)(p( P2 , S2 )\ J° c (0)|n(pi,3i)) 



The spin-averaged square of the amplitude Tbhvu can then be written as: 



|Tb#„„| 2 = 8waG F L V g H H^ , 



(147) 



(148) 



with the leptonic and hadronic tensors given by: 



L 



BH 



-Trj 

8 



(7" ^2 + 2^)7^(1-75) 

2(k'-q 2 ) 



7 /3 (1 - 75) (&7m + 2A£) 



2{k' -q 2 ) 



[k"q% + tflg - 9^ (*•«&)- if^KQ*, 



H ™ = \ E (p ^ S2 )i J - cc (°) i" (pi- s i)> (p (p 2 ' S2 )i J c (°) i n (pi. s i))* ' 



(149) 
(150) 



respectively. The matrix element of the weak charged transition current between the nucleon states is defined as: 



1 



(p (P2, «a)| #° (0) |n (pi, si)) = (p (p 2 , s 2 )| </> P (0) 7*5 ( X - 75 ) V>r» (0) |n (pi, 



(151) 



Using the isospin symmetry relation in Eq. I|130p between the flavor nondiagonal and diagonal nucleon matrix elements, 
we can parametrize the vector part of the matrix element in Eq. (j 1 5 1 j) in terms of the Dirac and Pauli form factors 
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for each quark flavor: 

(p(P2,s 2 )\4> p {0)j u ^n (0) \n(pi,si)) = (p(p 2 ,s 2 )\ 'tpu (0)luipu (0) \p(pi,si)) - (p(p2,S2)\ipd(0)j v ipd (0) \p(pi-,si)) 



u(p 2 ,s 2 ){ [Flu (t) - F ld (t)] lv - [F 2 „ (t) - F 2d (t)} 



2M 



u(pi,si) 



In addition, for the axial vector current part we have: 

(p(P2,S2)|^p(0)7,/75^T l (0)|n(pi,si)) = u(p 2 ,s 2 ) 9a (*) 7*75 ~9p (*) 



75 ?V 



2M 



u(pi,si) 



(152) 



(153) 



where the t dependence of the axial and pseudoscalar from factors follows the parametrizations in Eqs. I|123p and 



9A(t) = 9a 1 r 



.gp(i) = 9A{t) 



AM 2 

m 2 — t 



(154) 
(155) 



The hadronic tensor here is given by Eq. I|113p . with the form factors corresponding to the weak neutral transition 
current replaced by those describing the weak charged current interaction. The squared amplitude in Eq. l|148j) can 
then be written as: 



-8iraG%Ql 



(2k + t + l)x 2 B y 2 



x {x 2 B y 2 [2k + 1 + r (Ak + 1)] [[F lu (t) - F ld (t)] + [F 2u (t) - F 2d (t)]f 

+y [2x B (1 - r) - 4 + x B y [2 - x B (1 - r) + k [4 - 2x B (1 - 4 M )]]] 

t 



[Fiu (t) - Fid (t)]' 



4Af 



f2 



[F 2M (t)-*M (*)]' 



+y [2x B (1 - r) - 4 + zbj/ [2 - x B (1 - r) + /t [4 - 2.t b (1 + 4/i) + 4a?B-r]]] g 2 (*) 



a^y 2 (2k + t + 1) 



[g A (t)+g P (t)Y 



AM 2 



9 2 P (t) 



+Ax B y [1 + t - y {2k + 1) (1 + x b t)} [[F lu (t) - F ld (t)} + [F 2u (t) - F 2d (t)}} g A (t)} 



(156) 



Finally, the interference contribution between the Compton and Bethe-Heitler amplitudes has the following struc- 
ture: 



l un = -2iraG 2 F \ Tr <^ #' 7 „ (1 - 7s ) 



7 /3 (1-75)(^27m + 2^) 
2(k'-q 2 ) 



Tr< 



fc + M)7#£(f(i+M) 



[F lu (t) - F ld (t)\ 7/J + [F 2u (i) - F 2d (t)] ^f- - 5A (t) 7/57 5 ~ .9P (t) T "' •' 



+Tr ^(1-75) r 



2M 

(7" ^ 2 + 2^)7" (l-7s)' 



2M 



Tr\ (ji+M)(fP) (fa + M) 



2(k'-q 2 ) 

[F lu (t) - F ld (t)} 7 „ - [F 2u (t) - F 2d (t)] - g A (t) 7,75 + 9P (*) ^ 

in analogy with Eq. (|126[) . or more explicitly, after performing the contractions: 

1 



(157) 



-Z2-KaG 2 F Q\ 



1 



2k(2k + t + l)x B y 3 
-(2- x B ) (1 - y) [2 - y (2 - y)] + 2K [s - :r B (2 - y) 3 - y [11 - y (6 - y)] 
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[F lu (t) - F ld (t)] 3? (W 



4M- 



[F 2 » (t) - F 2d (t)] 3? (£ + c ) + g A (t) 3? 
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9 a (t) 
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9P (t) 



(158) 



Due to current conservation, the total amplitude, 

T^7i — ±Qi/ n -f" T BHvn i 



(159) 



is transverse with respect to the momentum q 2 of the outgoing real photon. In other words, by replacing the 
polarization vector e* (q 2 ) in Eqs. (( 1 2Qj) and (|147p by q 2u , we find for the Compton contribution: 



V2 



(Qu - Qd) |e| G F u (k') 7„ (1 - 75) u (k) u (p 2 , s 2 ) 



[F 2u (t) F 2d (t)] 17 J - g A (t) 7 ^ 75 + .gp (*) 



[Fiu (t) - F ld (t)\ 1 " 



AM 



2M 



u(pi,sx), 



(160) 



while for the Bethe-Heitler part we have: 



^BHv 



^ |e| G F u (k') Y (1 - 75) u (k) u { P2l s 2 ) 



[Fin (t) - Fi d (t)} l„ - [F 2u (t) - F 2d (t)} 



2M 



(161) 



Clearly, summing the expressions l|160p and l|16ip gives zero, so that the electromagnetic gauge invariance is explicitly 
satisfied. 
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E. Electromagnetic scattering 



For completeness (and for comparison with existing results [30]), we also consider the standard electromagnetic 
DVCS process on a proton target. The T-matrix for the pure Compton process is given by: 



1i 



(162) 



where 7™ is the reduced electromagnetic virtual Compton scattering amplitude, which can be easily reproduced 



from Eqs. I|39[) and (j43j) by discarding all terms accompanied with the "minus" OFPDs, and replacing the weak vector 
charge by the electric charge, c v — > Qf. One should also note the additional factor of 2 in the denominator of the 
amplitude due to the structure of the vertex qqZ°. The amplitude can then be written as: 
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( 92 • r) ; 
2Af 



where the corresponding integrals, known as Compton form factors, are given by: 

dx 



H+ (x, C,t) = Y, Q2 f J \ dx H f ( x > *) 
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(163) 

,(164) 
, (165) 
,(166) 
• (167) 



Here the spin- averaged square of the Compton T-matrix, Tc, can be written as: 



rC rrvp 



J4 ^uff^C > 



where the electron tensor, after neglecting the electron mass, becomes: 

L% = 2[k v k'p + k f} k' u -g V [ } (k-k')] 
and the electromagnetic hadronic tensor is: 

TTVp _ rr^lV ( \ * 

n c — 2 EM \ » EM ) 



(168) 



(169) 



-C\EM 



,"/8 



(P-Q2 



M 2 



{P ■ 92) 



1 - 



4M 2 



v 
<l2l2 



(170) 
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(171) 
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The contraction of tensors then yields: 



—2 _ (47ra) 

|Tcl " ~QW 2 



i + (i - 2/)' 



C\EM- 



(172) 



As in electron scattering via the ^°-exchange, the Bethe-Heitler contribution emerges from both diagrams (b) and 
(c) in Fig. G3 
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u(k)(p(p 2 , S2 )\J* M (0)\p( PuSl )), (173) 



where the proton matrix element of the electromagnetic transition current is parametrized in terms of the usual Dirac 
and Pauli proton form factors: 



(p( P2 ,s 2 )\J* M (0)\p( PllSl )) = u{p 2 ,s 2 ) 
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2M 



u(pi,si) 



Consequently one can write the Bethe-Heitler squared T-matrix as: 

(47TQ) jV p uBII 
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with the electron tensor: 
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(176) 



and the hadronic tensor: 
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The squared Bethe-Heitler amplitude is then found to be: 
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Finally, for the Compton-Bethe-Heitler interference part, we find: 
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It is worth noting that the analytic expressions for the Compton and interference contributions, see Eqs. I|172jl and 
(|179p . explicitly agree with the results of Ref. ^30] . As for the Bethe-Heitler contribution, one can recover the Bethe- 
Heitler result of Ref. [13] by substituting Eq. (|115p into the expression (|178|) . 

Having presented the formulas for cross sections for the electromagnetic and weak neutral and charged current 
DVCS, in the next section we compute these cross sections numerically, using the model for OFPDs in Sec. IIIIBl 



F. Cross section results 



The angular dependence of the Compton contribution to the unpolarized differential cross section for the weak 
neutral and charged DVCS processes is illustrated in Fig. [71 The results are presented for the region #b 7 < 12°, 
which corresponds to taking the invariant momentum transfer — t < — 1 GeV 2 (recall that the requirement for DVCS 
is that —t should be much smaller than Q\ = 2.5 GeV 2 ). The cross sections are observed to be of the same order in 
magnitude, and fall off smoothly with increasing 6*b 7 . For neutrino scattering, the vn charged current cross section 
is larger than the vp neutral current cross section. Furthermore, in the weak neutral current sector, the cross section 
is considerably larger for neutrinos rather than electrons. The latter reflects the difference in the structure of the 
leptonic tensors lu^l and L^l in Eqs. ([97| and <|104[1 . respectively, since the weak neutral hadronic tensor H^ N in 
Eq. i(98|) is the same in both cases. 
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Figure 7: Compton scattering cross sections for the weak DVCS processes corresponding to neutrino-proton scattering via the 
weak neutral current (solid), electron-proton scattering via the weak neutral current (dashed), and neutrino-neutron scattering 
via the weak charged current (dotted). The cross sections are plotted as a function of the scattering angle #b 7 between the 
incoming virtual weak boson and outgoing real photon in the target rest frame for Q? = 2.5 GeV 2 and xb = 0.35 with an 
ui = 20 GeV lepton beam. 

To investigate the magnitude of the Compton contributions relative to the corresponding Bethe-Heitler backgrounds, 
we plot both cross sections together on a logarithmic scale in Figs. [8] and El Note that for the Bethe-Heitler cross 
section, one should expect poles when the scattering angle 6*b 7 coincides with 4> — 20.2° and <j>' = 25.3°, which 
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corresponds to the outgoing real photon being collinear with either the incoming or scattered lepton. In both cases, 
the Compton cross section is significantly larger than the Bethe-Heitler background, especially at small angles. This 
is to be contrasted with the electromagnetic current case, Fig. [TOl in which the Bethe-Heitler contribution is larger 
than the Compton for 6* 77 > 3°. 
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Figure 8: Compton (dashed), Bethe-Heitler (dotted), magnitude of interference (dash-dotted) and total (solid) cross sections 
for electron-proton scattering via the weak neutral current plotted as a function of the scattering angle between the incoming 
virtual weak boson and outgoing real photon in the target rest frame for Q\ — 2.5 GeV 2 and xb = 0.35, with an it) = 20 GeV 
electron beam. 
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Figure 9: As in Fig. [8] but for the neutrino-neutron charged current cross section. 
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Figure 10: As in Fig. [9l but for the electromagnetic electron-proton cross section. 

The relative importance of the various contributions is even more graphically illustrated in Fig. Qj], where we plot 
the ratio of the Compton and Bethe-Heitler cross sections for the weak neutral current (solid) , weak charged current 
(dashed), and electromagnetic (dotted) DVCS. In the forward direction (e.g. for 6*s 7 < 4° with the ratio —t/Q\ < 0.1), 
the Bethe-Heitler contribution is strongly suppressed by a factor more than 100 compared to the Compton contribution 
for both weak neutral and weak charged current scattering. This is in contrast with the electromagnetic current case, 
where the ratio between the cross sections remains of order unity. Although based on a rather simple model of OFPDs, 
these results should be helpful in providing some guidance for future neutrino scattering experiments. 
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Figure 11: Ratio of Compton and Bethe-Heitler cross sections for weak neutral (solid), weak charged (dashed) and electromag- 
netic (dotted) DVCS, as a function of the scattering angle between the incoming boson B = Z°, W + or 7 and the outgoing 
real photon in the target rest frame, for Q\ = 2.5 GeV 2 , xb — 0.35 and lo = 20 GeV. 



IV. CONCLUSIONS 



In this work we have presented a comprehensive analysis of DVCS-like processes induced by the exchange of weak 
vector bosons. Weak DVCS is an important new tool for studying the quark structure of the nucleon, providing 
complementary information to that available through electromagnetic probes. 
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As in elastic form factor studies, or in deep inelastic scattering, weak DVCS provides access to different combinations 
of GPDs to those which can appear in ordinary DVCS. The flavor dependence of the weak charges, for example, 
means that one has more sensitivity to d quarks in the proton than in electromagnetic scattering. Furthermore, the 
V — A nature of the weak interaction allows both C-odd and C-even combinations of GPDs to be measured, thereby 
providing a direct separation of the valence and sea content of GPDs. Since DVCS is sensitive to spin-averaged 
as well as spin-dependent GPDs, this feature means that one can extract spin-dependent valence and sea quark 
distributions without polarizing the target. The only other means of obtaining spin-dependent C-odd distributions 
would be through neutrino-polarized nucleon scattering, which is prohibitive, however, because of the large quantities 
of nuclear material which would need to polarized. An additional feature of weak DVCS is that it enables one to 
access GPDs that are nondiagonal in quark flavor, such as those associated with the neutron-to-proton transition. 

In the present study, we have derived the weak virtual Compton scattering amplitude for the neutral current, in 
both neutrino and charged lepton scattering, as well as for the charged current. The amplitudes have been calculated 
in the leading, twist-2 approximation using the light-cone expansion of the current product, in terms of QCD string 
operators in coordinate space. 

To quantify our results, and to explore the feasibility of measuring weak DVCS cross sections, we have used a simple, 
factorized model to estimate the cross sections in kinematics relevant to future high-intensity neutrino experiments. 
In contrast to the standard electromagnetic DVCS process, we find that at small scattering angles the Compton 
signal is enhanced relative to the corresponding Bethe-Heitler contribution. This should make contamination from 
the Bethe-Heitler backgrounds less of a problem when extracting the weak DVCS signal. 

While the current model analysis has been exploratory, in future one can use more elaborate models for nucleon 
GPDs, including sea quark effects, and contributions from the plus and minus distributions separately. In the small- 
\t\ region of DVCS kinematics, it would also be of interest to examine further the contribution from the pion pole, 
through the Ef distribution. Finally, it will be necessary to extend the approach to include twist-3 terms in order to 
apply the formalism at moderate energies, where the suppression of higher twist contributions is not guaranteed. The 
results presented here provide an important starting point in realizing the program of extracting GPDs in neutrino 
scattering. 
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